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CN , The masses of the excited heavy tetraquarks with hidden charm are calculated 

bJQI within the relativistic diquark-antidiquark picture. The dynamics of the light quark 
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< 



X 



in a heavy-light diquark is treated completely relativistically. The diquark structure 
is taken into account by calculating the diquark-gluon form factor. New experimental 
2^ . data on charmonium-like states above open charm threshold are discussed. The 

obtained results indicate that X(3872), y(4260), ^(4360), Z(4248), Z(4433) and 
' f— I 'I y(4660) could be tetraquark states with hidden charm. 

ji^! PACS numbers: 12.40.Yx, 14.40.Gx, 12.39.Ki 

Recently, significant experimental progress has been achieved in charmonium spec- 
troscopy. Several new states, such as X(3872), F(4260), r(4360), r(4660), Z(4248), 
' Z(4430), etc., were observed [l| which cannot be simply accommodated in the quark- 
^ ■ antiquark (cc) picture. These states and especially the charged ones can be considered 
as indications of the possible existence of exotic multiquark states [2I, Isj. In our papers 
\ 0, H] we calculated masses of the ground state heavy tetraquarks in the framework of the 

00 ■ relativistic quark model based on the quasipotential approach in quantum chromodynamics. 

Here we extend this analysis to the consideration of the excited tetraquark states with hid- 
! den charm. As previously, we use the diquark-antidiquark picture to reduce a complicated 
relativistic four-body problem to the subsequent two more simple two-body problems. The 
first step consists in the calculation of the masses, wave functions and form factors of the 
diquarks, composed from light and heavy quarks. At the second step, a heavy tetraquark is 
5^ \ considered to be a bound diquark-antidiquark system. It is important to emphasize that we 
do not consider the diquark as a point particle but explicitly take into account its structure 
by calculating the form factor of the diquark-gluon interaction in terms of the diquark wave 
functions. 

In the quasipotential approach and diquark-antidiquark picture of heavy tetraquarks 
the interaction of two quarks in a diquark and the diquark-antidiquark interaction in a 
tetraquark are described by the diquark wave function (\E'rf) of the bound quark-quark state 
and by the tetraquark wave function {^t) of the bound diquark-antidiquark state, respec- 
tively. These wave functions satisfy the quasipotential equation of the Schrodinger type 
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where the relativistic reduced mass is 
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Here, M = Ei + E2 the bound-state mass (diquark or tetraquark), mi 2 are the masses 
of quarks (g and Q) which form the diquark or of the diquark {d) and antiquark {d') which 
form the heavy tetraquark (T), and p is their relative momentum. In the center-of-mass 
system the relative momentum squared on mass shell reads 
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The kernel l^(p, q; M) in Eq. ([T]) is the quasipotential operator of the quark-quark or 
diquark-antidiquark interaction. It is constructed with the help of the off-mass-shell scatter- 
ing amplitude, projected onto the positive-energy states. In the following analysis we closely 
follow the similar construction of the quark-antiquark interaction in mesons which were ex- 
tensively studied in our relativistic quark model |1, 0] • For the quark-quark interaction in a 
diquark we use the relation Vqq = Vgg/2 arising under the assumption of an octet structure of 
the interaction from the difference in the qq and qq colour states. An important role in this 
construction is played by the Lorentz structure of the confining interaction. In our analysis 
of mesons, while constructing the quasipotential of the quark-antiquark interaction, we as- 
sumed that the effective interaction is the sum of the usual one-gluon exchange term and a 
mixture of long-range vector and scalar linear confining potentials, where the vector confin- 
ing potential contains the Pauli term. We use the same conventions for the construction of 
the quark-quark and diquark-antidiquark interactions in the tetraquark. The quasipotential 
is then defined as follows [Z, iSi] . 

(a) For the quark-quark (Qq) interactions, V{p, q; M) reads 



V{p, q; M) = Ui{p)u2{-p)V{p, q; M)ui{q)u2{-q), 
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with 



V(p,q;M) = - 



Here, is the QCD coupling constant; D^i, is the gluon propagator in the Coulomb gauge, 
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and k = p — q; 7^ and u{p) are the Dirac matrices and spinors. 
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with e{p) = + m^. 
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The effective long-range vector vertex of the quark is defined 0] by 

~k = {0,k), 



IK 
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where k is the Pauh interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In configuration space the vector and scalar confining potentials in the 
nonrelativistic limit 191 reduce to 



with 



Konf (r) = Konf (0 + Konf (0 = Av + B 
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where e is the mixing coefficient. 

(b) For the diquark-antidiquark {dd') interaction, V^(p, q; M) is given by 

+rd{p)rd'{n \Jd;,Ji;>vLi{^) + Konf(k)l mq)MQ') 



where {d{P)\J ^\d{Q)) is the vertex of the diquark-gluon interaction which takes into account 
the finite size of the diquark and is discussed below P*^') = (-E^o, ±p) and Q^'^ = (-E^o, 
Ed = (M^ - M| + M|)/(2M) and E^, = (M^ - Mj + M|)/(2M) " . 

The diquark state in the confining part of the diquark-antidiquark quasipotential (fTTi) is 
described by the wave functions 



1 

£d{p) 



for a scalar diquark, 

for an axial-vector diquark. 



where the four-vector 
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is the polarization vector of the axial- vector diquark with momentum p, Edip) = -y/P 
and £d(0) = (0, £d) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 



J, 



2jEdEd 
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Here, the antisymmetric tensor EJ;^ is defined by 



for a scalar diquark, 

for an axial-vector diquark. 
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TABLE I: Masses M and form factor parameters of charmed diquarks. S and A denote scalar and 
axial vector diquarks which are antisymmetric [• • •] and symmetric {• • •} in flavour, respectively. 



Quark 
content 


Diquark 
type 


M 
(MeV) 


(GeV) 


c 

(GeV2) 


[c, q] 


S 


1973 


2.55 


0.63 


{c, q} 


A 


2036 


2.51 


0.45 


[c, s] 


S 


2091 


2.15 


1.05 


{c,s} 


A 


2158 


2.12 


0.99 



and the axial-vector diquark spin Sd is given by {Sd;k)ii = —i^ku] fJ'd is the total chromo- 
magnetic moment of the axial-vector diquark. 

The constituent quark masses rric = 1.55 GeV, = = 0.33 GeV, = 0.5 GeV and 
the parameters of the linear potential A = 0.18 GeV^ and B = —0.3 GeV have values typical 
in quark models. The value of the mixing coefficient of vector and scalar confining potentials 
e = —1 has been determined from the consideration of charmonium radiative decays [6| and 
the heavy-quark expansion p^]. The universal Pauli interaction constant k = —1 has been 
fixed from the analysis of the fine splitting of heavy quarkonia ^Pj - states [oj. In this 
case, the long-range chromomagnetic interaction of quarks vanishes in accordance with the 
fiux-tube model. 

At the first step, we calculate the masses and form factors of the heavy-light diquark. 
As it is well known, the light quarks are highly relativistic, which makes the v/c expan- 
sion inapplicable and thus, a completely relativistic treatment of the light quark dynamics 
is required. To achieve this goal, we closely follow our consideration of diquarks in heavy 
baryons and adopt the same procedure to make the relativistic potential local by replacing 

ei,2(p) = \J^\,2 + ~^ -^1,2 = — "^2,1 + Ta\ 2)l'^M- Solving numerically the quasipo- 
tential equation ([T]) with the complete relativistic potential, which depends on the diquark 



mass in a complicated highly nonlinear way [Uj , we get the diquark masses and wave func 



tions. In order to determine the diquark interaction with the gluon field, which takes into 
account the diquark structure, we calculate the corresponding matrix element of the quark 
current between diquark states. Such calculation leads to the emergence of the form factor 
F{r) entering the vertex of the diquark-gluon interaction This form factor is expressed 
through the overlap integral of the diquark wave functions. Our estimates show that this 
form factor can be approximated with a high accuracy by the expression 



(16) 



The values of the masses and parameters ^ and C, for heavy-light scalar diquark [■ ■ ■] and 
axial vector diquark {■ ■ ■} ground states are given in Table [H 

At the second step, we calculate the masses of heavy tetraquarks considered as the bound 
states of a heavy-light diquark and antidiquark. For the potential of the diquark-antidiquark 
interaction fill I) we get |5| 
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where 



y 4 Fi(r)F2(r) 

Vcoui(r) = --q;^ 

6 r 



is the Coulomb-hke one-gluon exchange potential which takes into account the finite sizes 
of the diquark and antidiquark through corresponding form factors Fi^2{t)- Here, 81,2 and 

L are the spin operators of diquark and antidiquark and the operator of the relative orbital 
angular momentum. In the following we choose the total chromomagnctic moment of the 
axial- vector diquark /x^ = 0. Such a choice appears to be natural, since the long-range 
chromomagnctic interaction of diquarks proportional to //^ then also vanishes in accordance 
with the flux-tube model. 

In the diquark-antidiquark picture of heavy tetraquarks both scalar S (antisymmetric in 
flavour {Qq)s=Q = [Qq]) and axial vector A (symmetric in flavour (Qq)s=i = {Qq}) diquarks 
are considered. Therefore we get the following structure of the {Qq){Qq') ground {IS) states 
(C is defined only for q — q'): 



Two states with J^^ = 0++: 



X(0++) = {Qq)s=o{Qq)s=^ 
X(0++') = {Qq)s=i{Q<i)s=i 



• Three states with J = 1: 



X(l++) = -j^[{Qq)s=i{Q(i)s=o + (Qg)5=o(Qg')5=i] 
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TABLE II: Masses of charm diquark-antidiquark ground {IS) states (in MeV) calculated in [J]. 5 
and A denote scalar and axial vector diquarks. 



State 


Diquark 




Mass 






content 


cqcq 


cscs 


cqcs 


0++ 


SS 


3812 


4051 


3922 


l+± 


{SA±SA)/V2 


3871 


4113 


3982 


0++ 


AA 


3852 


4110 


3967 


1+- 


AA 


3890 


4143 


4004 


2++ 


AA 


3968 


4209 


4080 



TABLE III: Thresholds for open charm decays and nearby hidden-charm thresholds. 
Channel Threshold (MeV) Channel Threshold (MeV) Channel Threshold (MeV) 





3729.4 




3936.2 




3832.9 


D+D- 


3738.8 


ry'J/V' 


4054.7 


D^Df 


3837.7 




3871.3 


DfDt^ 


4080.0 




3975.0 




3872.7 




4116.4 




3976.7 




3879.5 


Dl+D*- 


4223.8 


K*^J/ij 


3988.6 


uiJ/tp 


3879.6 






K*^J/iIj 


3993.0 




4013.6 








4118.8 



^(1+") = ^[ms=oiQq')s=i - {Qq)s=iiQq)s=o] 
X(l+^') = {Qq)s=i{Qq')s=i 

• One state with J^"^ = 2++: 

X(2++) = {Qq)s=i{Qq')s=i. 

The orbitally excited (IP, ID...) states are constructed analogously. As we find, a very 
rich spectrum of tetraquarks emerges. However the number of states in the considered 
diquark-antidiquark picture is significantly less than in the genuine four-quark approach. 

The diquark-antidiquark model of heavy tetraquarks predicts the existence of a flavour 
SU{3) nonet of states with hidden charm or beauty {Q = c,b): four tetraquarks [{Qq){Qq), 
q = u,d] with neither open or hidden strangeness, which have electric charges or ±1 and 
isospin or 1; four tetraquarks [{Qs){Qq) and {Qq){Qs), q = u,d] with open strangeness 
{S = ±1), which have electric charges or ±1 and isospin |; one tetraquark {Qs){Qs) 
with hidden strangeness and zero electric charge. Since we neglect in our model the mass 
difference of u and d quarks and electromagnetic interactions, the corresponding tetraquarks 
will be degenerate in mass. A more detailed analysis [12] predicts that the tetraquark mass 
differences can be of a few MeV so that the isospin invariance is broken for the {Qq){Qq) 
mass eigenstates and thus in their strong decays. The (non)observation of such states will 
be a crucial test of the tetraquark model. 
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TABLE IV: Masses of charm diquark-antidiquark excited IP, 25" states (in MeV). S and A denote 
scalar and axial vector diquarks; S is the total spin of the diquark and antidiquark. (C is defined 



only for q = q'). 


State 


Diquark 






Mass 






content 


S 


cqcq 


cscs 


cqcs 


IP 












1— 


SS 





4244 


4466 


4350 


o-± 


{SA±SA)/V2 


1 


4269 


4499 


4381 


l-± 


{SA±SA)/V2 


1 


4284 


4514 


4396 




{SA±SA)/V2 


1 


4315 


4543 


4426 


1— 


AA 





4350 


4582 


4461 


0-+ 


AA 


1 


4304 


4540 


4419 


1-+ 


AA 


1 


4345 


4578 


4458 


2-+ 


AA 


1 


4367 


4598 


4478 


1 — 


AA 


2 


4277 


4515 


4393 


2— 


AA 


2 


4379 


4610 


4490 


3— 


AA 


2 


4381 


4612 


4492 


2S 












0++ 


SS 





4375 


4604 


4481 


l+± 


{SA±SA)/V2 


1 


4431 


4665 


4542 


0++ 


AA 





4434 


4680 


4547 


1+- 


AA 


1 


4461 


4703 


4572 


2++ 


AA 


2 


4515 


4748 


4625 



The calculated masses of the heavy tetraquark ground (IS*) states and the corresponding 
open charm thresholds are shown in Tables [ITl HHi Note that most of the tetraquark 
states were predicted to lie either above or only slightly below corresponding open charm 
thresholds. In Table \IV\ IVl we give our predictions for the orbitally and radially excited 
tetraquark states with hidden charm. Excitations only of the diquark-antidiquark system 
are considered. A very rich spectrum of excited tetraquark states is obtained. 

In Table I VI I we compare our results (EFG) for the masses of the ground and excited 



charm diquark-antidiquark bound states with the predictions of Refs. [12|, ll3|, Il 4 Il5 | and 



with the masses of the rece ntly observed highly-excited charmonium-like states 



18, m 



21, 22 



25 



We assume that the excitations occur only between the 
bound diquark and antidiquark. Possible excitations of diquarks are not considered. Our 
calculation of the heavy baryon masses supports such a scheme jllli]. In this table we give 
our predictions only for some of the masses of the orbitally and radially excited states 
for which possible experimental candidates are observed. The differences in some of the 
presented theoretical mass values can be attributed to the substantial distinctions in the 
used approaches. We describe the diquarks dynamically as quark-quark bound systems and 
calculate their masses and form factors, while in Refs. (121. 13 . III . 15] they are treated only 
phenomenologically. Then we consider the tetraquark as purely the diquark-antidiquark 
bound system. In distinction, Maini et al. consider a hyperfine interaction between all 
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TABLE V: Masses of charm diquark-antidiquark excited ID, 2P states (in MeV). S and A denote 
scalar and axial vector diquarks; S is the total spin of the diquark and antidiquark. 



State 


Diquark 






Mass 






content 


s 


cgcg 


cscs 


cqcs 


ID 












2++ 


SS 





4506 


4728 


4611 


l+± 


{SA±SA)/V2 


1 


4553 


4779 


4663 


2+± 


{SA±SA)/y/2 


1 


4559 


4785 


4670 


3+± 


{SA±SA)/V2 


1 


4570 


4794 


4680 


2++ 


AA 





4617 


4847 


4727 


1+- 


AA 


1 


4604 


4835 


4714 


2+- 


AA 


1 


4616 


4846 


4726 


3+- 


aa 


1 


4624 


4852 


4733 


0++ 


aa 


2 


4582 


4814 


4692 


1++ 


AA 


2 


4593 


4825 


4703 


2++ 


AA 


2 


4610 


4841 


4720 


3++ 


AA 


2 


4627 


4855 


4736 


4++ 


AA 


2 


4628 


4856 


4738 


2P 












1— 


SS 





4666 


4884 


4767 


o-± 


{SA±SA)/^/2 


1 


4684 


4909 


4792 


l-± 


{SA±SA)/^/2 


1 


4702 


4926 


4810 


2-± 


{SA±SA)/V2 


1 


4738 


4960 


4845 


1— 


AA 





4765 


4991 


4872 


0-+ 


AA 


1 


4715 


4946 


4826 


1-+ 


AA 


1 


4760 


4987 


4867 


2-+ 


AA 


1 


4786 


5011 


4892 


1— 


AA 


2 


4687 


4920 


4799 


2— 


AA 


2 


4797 


5022 


4903 


3— 


AA 


2 


4804 


5030 


4910 



quarks which, e.g., causes the spHtting of l"*"^ and 1'^^ states arising from the SA diquark- 
antidiquark compositions. From Table IVII we see that our dynamical calculation supports 
the assumption [l2|] that X(3872) can be the axial vector tetraquark state composed 
from the scalar and axial vector diquark and antidiquark in the relative IS" state. Recent 
Belle and BaBar results indicate the existence of a second X(3875) particle a few MeV above 
X(3872). This state could be naturally identified with the second neutral particle predicted 



by the tetraquark model 13|. On the other hand, in our model the lightest scalar 0^"*" 
tetraquark is predicted to be above the open charm threshold DD and thus to be broad, 
while in the model fl3\ it lies a few MeV below this threshold, and thus is predicted to 



be narrow. Our 2"*"+ tetraquark also lies higher than the one in Ref.[12|, thus making the 



interpretation of this state as F(3943) less probable, especially if one averages the original 
Belle result with the recent BaBar value which is somewhat lower. 
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TABLE VI: Comparison of theoretical predictions for the masses of the ground and excited charm 
diquark-antidiquark states (in MeV) and possible experimental candidates. 



State 


Diquark 




Theory 




Experiment 




content 


EFG 


[12,13,14] 


[15] {cscs) 


state 


mass 


IS 














0++ 


SB 


3812 


3723 








1 ++ 


[bA + bA)ly 2, 


3871 


3872' 


1 

1 


' X(3872) f 
X(3876) \ 


3871.4 ±0.6 [1] 
3875.2 ± 0.7+?i [1] 


1+- 


{SA - SA)/^/2 


3871 


3754 








0++ 


AA 


3852 


3832 








1+- 


AA 


3890 


3882 








2++ 


AA 


3968 


3952 




y(3943) 1 


3943 ±11 ±13 [16] 
3914 \17] 


IP 














I— 


SS 


4244 




4330±70 


y(4260) 


f 4259 ± 8j:g [18] 
1 4247 ± fl9l 


1- 


SS 

{SA ± SA)/^/2 


4244 ' 
4267 ^ 


1 

f 




Z(4248) 


A2A8ttttll' [20] 


1~ 
1~ 


{SA-SA)/V2 
AA 


4284 ' 
4277 ^ 


1 

\ 




y(4260) 


4284t}^±4 [21] 


1— 


AA 


4350 






I y'ioytyj ) 


f 4361 ± 9 ± 9 [22] 
\ 4324 ± 24 [23] 


25 














1+ 
0+ 


{SA±SA)/y/2 
AA 


4431 ' 
4434 


\ 




Z(4430) 


4433±4±2 [24] 


1+ 


AA 


4461 


~ 4470 








2P 














1— 


SS 


4666 




1 


' y(4660) 
_ X(4630) 


' 4664 ± 11 ± 5 [22] 
^ 46341^1^ [25] 



t \ 



input 



The recent discovery in the initial state radiation at 5-factories of the y(4260), y(4360) 
and y (4660) indicates an overpopulation of the expected charmonium 1 states llSl. Il9l. 



21, 22, 23 



Maini et al. [15|] argue that y(4260) is the 1 IP state of the charm-strange 
diquark-antidiquark tetraquark. We find that y(4260) cannot be interpreted in this way, 
since the mass of such ([cs]5'=o[cs]5=o) tetraquark is found to be ~ 200 MeV higher. A more 
natural tetraquark interpretation could be the 1 IP state ([cg]5=o[cg]5'=o) (SS) which 
mass is predicted in our model to be close to the mass of y(4260) (see Table IVll) . Then 
the y(4260) would decay dominantly into DD pairs. The other possible interpretations 
of y(4260) are the 1 IP states of {SA — SA)/\/2 and AA tetraquarks which predicted 
masses have close values. These additional tetraquark states could be responsible for the 
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mass difference of y(4260) observed in different decay cliannels. As we see from Table IVlt 
tlie recently discovered resonances F (4360) and y (4660) in the e'^e~ 7r^n~ip' cross section 
can be interpreted as the excited 1 IP (AA) and 2P (SS) tetraquark states, respectively. 
The peak X(4630) very recently observed by Belle in e~^e^ ^t-^c consistent with 

a 1 resonance y(4660) and therefore has the same interpretation in our model. 

Recently the Belle Collaboration reported the observation of a relatively narrow enhance- 
ment in the Tf^ip' invariant mass distribution in the B — >■ Kn^ip' decay [H, 24|. This new 
resonance, Z"''(4430), is unique among other exotic meson candidates, since it is the first 
state which has a non-zero electric charge. Different theoretical interpretations were sug- 
gested [l|. Maiani et al. [3] give qualitative arguments that the Z+(4430) could be the 
first radial excitation (25) of a diquark-antidiquark X^(l^ ; 15) state {AA) with mass 
3882 MeV. Our calculations indicate that the Z+(4430) can indeed be the 1+ 25" [cM][c(i] 
tetraquark state. It could be the first radial excitation of the ground state (SA — SA)/^/2, 
which has the same mass as X(3872). The other possible interpretation is the O"*" 2S [cu][cd] 
tetraquark state (AA) which has a very close mass. Measurement of the Z"'"(4430) spin will 
discriminate between these possibilities. 

Encouraged by this discovery, the Belle Collaboration performed a study of — >■ 
K~n~^-Xci and observed a double peaked structure in the n~^Xci invariant mass distribu- 
tion [20]. These two charged hidden charm peaks, Z(4051) and Z(4248), are explicitly 
exotic. We find no tetraquark candidates for the former, Z(4051), structure. On the other 
hand, we see from Table IVl] that Z(4248) can be interpreted in our model as the charged 
partner of the 1~ IP state SS or as the 0~ IP state of the {SA ± SA)/\/2 tetraquark. 

In summary, we calculated the masses of excited heavy tetraquarks with hidden charm 
in the diquark-antidiquark picture. In contrast to previous phenomenological treatments, 
we used the dynamical approach based on the relativistic quark model. Both diquark and 
tetraquark masses were obtained by numerical solution of the quasipotential wave equation 
with the corresponding relativistic potentials. The diquark structure was taken into account 
in terms of diquark wave functions. It is important to emphasize that, in our analysis, we did 
not introduce any free adjustable parameters but used their values fixed from our previous 
considerations of heavy and light hadron properties. It was found that the X(3872), Z(4248), 
F(4260), y(4360), Z(4430) and F(4660) exotic meson candidates can be tetraquark states 
with hidden charm. 
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